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Abstract -Critical phenomena can show unusual phase diagrams when defined in complex net- 
work topologies. The case of classical phase transitions such as the classical Ising model and the 
percolation transition has been studied extensively in the last decade. Here we show that the 
phase diagram of the Bose-Hubbard model, an exclusively quantum mechanical phase transition, 
also changes significantly when defined on random scale-free networks. We present a mean-field 
calculation of the model in annealed networks and we show that when the second moment of 
the average degree diverges the Mott-insulator phase disappears in the thermodynamic limit. 
Moreover we study the model on quenched networks and we show that the Mott-insulator phase 
disappears in the thermodynamic limit as long as the maximal eigenvalue of the adjacency matrix 
diverges. Finally we study the phase diagram of the model on Apollonian scale-free networks that 
can be embedded in 2 dimensions showing the extension of the results also to this case. 



Introduction. — Recently great attention [Ud] has 
been addressed to critical phenomena unfolding on com- 
plex networks. In this context it has been observed that 
the topology of the networks might significantly change 
the phase diagram of dynamical processes. For exam- 
ple when networks have a scale-free degree distribution 
P(k) ~ fc~ A and the second moment (k 2 ) diverges with the 
network size, i.e. A S (2, 3] the Ising model (3]-[6], the per- 
colation phase transition [TIE] and the epidemic spreading 
dynamics on annealed networks [5] are strongly affected. 
Moreover the spectral properties of the networks drive the 
epidemic spreading on quenched networks [TCTIfTT) , the syn- 
chronization stability |12ll3j . the critical behavior otO(N) 
models [HI [15] and the critical fluctuations of an Ising 
model on spatial scale- free networks |16| . Quantum criti- 
cal phenomena also might depend on the topology of the 
underlying lattice as it has been shown for Bosc-Einstein 
condensation in heterogeneous networks |17j . Although 
large attention has been devoted to classical critical phe- 
nomena on scale-free networks, the behavior of quantum 
critical phenomena on scale-free networks has just started 
to be investigated. In particular the Anderson local- 
ization |18H19| was studied in complex networks show- 
ing that by modulating the clustering coefficient of the 



network one might induce localization transition in scale- 
free networks. Moreover, attention has been addressed to 
quantum processes on Apollonian networks 20,21 , which 
provide an example of scale-free networks embedded in 
two dimensions. The quantum processes investigated are 
the Hubbard model [35], the free electron gas within the 
tight-binding model [53], and the topology induced Bosc- 
Einstein condensation in Apollonian networks [24]. The 
study of quantum phase transitions on these networks 
has attracted attention particularly in recent years moti- 
vated by the creation of a new self-similar macromolecule 
- a nanometer-scale Sierpinski hexagonal gasket |25j . Re- 
cently [26] it has been shown that the Random Transverse 
Ising model is strongly affected by a scale-free network 
topology of the underlying networks on which it is defined 
and in particular by the second moment of the degree dis- 
tribution (k 2 ). Indeed the critical temperature for the 
onset of the disordered phase is infinite if this second mo- 
ment diverges and the network is scale-free with power-law 
exponent A < 3. 

In this paper we investigate a critical process with no 
classical equivalent, the Bose-Hubbard model on complex 
networks [271430] . In the framework of ultracold atom 
physics nowadays defect free potentials have been con- 
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structed and the phenomenology of systems described by 
the Bose-Hubbard model experimentally reproduced [3"T] . 
By superimposing different lattices with incommensurate 
lattice constants, disordered systems of ultracold atoms 
have been recently experimentally investigated [32]. As 
the experimental technology advances it might become 
possible to investigate the role of topological network 
complexity in the phase diagram of the Bose-Hubbard 
model. Moreover this complex topology might be effec- 
tively present in complex granular materials and might 
affect insulator-superconductor phase transitions [33] in 
these systems therefore a full account of the consequence 
of complex topologies might turn out to shed some light in 
the phase diagram of these complex materials. In fact the 
model defined on complex networks might provide a use- 
ful mean-field approximation to real disordered granular 
materials that captures essential features of their hetero- 
geneity. Finally the Hubbard model is a theoretical model 
that has applications far beyond condensed matter physics 
(SUES] and investigating its properties on scale-free net- 
works might stimulate further applications to other fields. 
Here we investigate this model when it is defined on a 
scale-free network topology. Recently different approaches 
have been suggested for the theoretical study of the Bose- 
Hubbard model. Here we cite field theoretic approxima- 
tions [28I36H38] . mean-field approximations, [29,30,39,40 
quantum Monte Carlo simulations J41H46] . quantum cav- 
ity methods |47j among many others. 

In this paper we characterize the phase diagram of the 
Bose-Hubbard model by mean-field approximation on an- 
nealed and quenched networks. The Bose-Hubbard model 
is described by the Hamiltonian 



H 



E 



u 



n,-(rii — 1) - i-m, 



'E' 



(i) 



i -j 



where the indices i,j = 1,...N indicate the nodes of the 
network. The network has an adjacency matrix r such 
that Tij = 1 if there is a link between node i and node 
j, otherwise ry = 0. The operator a\(ai) creates (annihi- 
lates) a boson at site i and n, = a\ai counts the bosons at 
site i. The parameter U represents the repulsive boson- 
boson interaction, t represents the hopping amplitude be- 
tween neighboring nodes while \x indicates the chemical 
potential. The hallmark of this Hamiltonian is a quantum 
phase transition between the Mott insulating phase and 
the supcrfluid phase originating from the competition be- 
tween the kinetic and the repulsive terms of the Hamilto- 
nian. We have to mention that the Mott insulating phase 
with vanishing compressibility is present strictly speaking 
only at zero temperature. At finite temperature thermal 
fluctuations induce a phase transition between the super- 
fluid and normal phase. 

In this paper we show by mean-field approximations 
that the phase diagram of the model defined on an an- 
nealed network depends on the second moment of the de- 
gree distribution (k 2 ). In particular, by the mean- field 



approximation, we found both for annealed and quenched 
networks that for scale-free networks with A < 3 the Mott 
insulator phase reduces with increasing network size, dis- 
appearing in the thermodynamic limit. Moreover we ob- 
serve differences between the model defined on a quenched 
network and an annealed network. In fact it is sufficient 
for a quenched random network to have diverging maxi- 
mum degree in order to reduce the Mott-insulator phase 
to zero in the thermodynamic limit. This demonstrates 
that complex networks might strongly perturb the phase 
diagram of quantum phase transitions. 

The paper is structured as follows: In section II we give 
the solution of the Bose-Hubbard model in annealed com- 
plex networks within the mean-field approximation. In 
section III we study numerically the phase diagram of the 
Bose-Hubbard model on quenched scale-free networks in 
the mean-field approximation and also on Apollonian net- 
works. Finally we give the concluding remarks. 

Mean-field solution of the Bose-Hubbard Model 
on annealed complex networks. — An annealed net- 
work evolves dynamically on the same time scale as the 
dynamical process occurring on it. During this process, 
links are created and annihilated but the expected degree 
of each node remains the same. Solving dynamical mod- 
els in annealed networks is usually straightforward. We 
are nevertheless in general not guaranteed that the phase 
diagram of the dynamical model on annealed scale-free 
networks will capture the essence of the dynamical model 
on quenched networks. Our strategy here will be first to 
study the Bose-Hubbard model in the annealed approxi- 
mation and then to study the model on quenched networks 
to validate the main conclusions of the paper. 

We consider the ensemble of uncorrelated networks in 
which we assign to each node a hidden variable 0$ from a 
distribution p(9) indicating the expected number of neigh- 
bors of a node. We consider the ensemble of networks for 
which the probability to draw a link between node i and 
j is given by p tj 



Pij = P(nj = l) 



(9)N- 



(2) 



In this ensemble the degree fcj of a node i is a Poisson 
random variable with expected degree hi — 6i. Therefore 
we will have 



W = w 

l$ 2 ) = (fc(fc-i)}. 



(3) 



where (. . .) indicates the average over the N nodes of the 
network and the overline in Eq. ([3]) indicates the average 
over the ensemble of the networks. We assume that the 
expected degree distribution of the network ensemble is 
given by 



V 



(9)=JV9- x e- e/ t 



(1) 



where A/" is a normalization constant and £ is an exponen- 
tial cut-off in the expected degree distribution. 
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In order to study the Bose-Hubbard model on annealed 
complex networks we consider the fully connected Hamil- 
tonian given by 



H = ^2—m(m 



1) - nm - tS^ pijCLia] 



(5) 



• ■j 



where in order to account for the dynamical nature of the 
annealed graph we have substituted the adjacency matrix 
element r^- in H given by Eq. ([T]) with the matrix clement 
Pij given by Eq. O Moreover we perform the mean-field 
approximation to the Bose-Hubbard model introduced in 
[39] by taking 

{Ha'j ~ (ai)aj + ai(a!j) - (ai)(a]} 

~ ipitt] + cui/jj - ipiipj (6) 

where V»« = (o») = (oj) ■ The Hamiltonian is then decom- 
posed in single site terms 



H = J2H t + {6)NtY 



(T) 



with Hi given by 

Hi = —rii(rii - 1) - pn t - t9 t j(a t + a\) (8) 

and with 7 indicating the order parameter of the supcrfluid 
phase, defined as 



7= P^ 9 * 



(9) 



In this mean-field picture the Hamiltonian therefore de- 
couples in single site (node) Hamiltonians Hi depending 
on the mean field order parameter 7. We can therefore 
write the single site (node) Hamiltonian as an unperturbed 
Hamiltonian plus an interaction depending on the param- 
eter 7, i.e. 



Hi = H 



(0) 



lOiVi 



(10) 



with 



H 



(0) 



U 



m(ni-l)- pn t 



Vi = t(<n + al) 



(2) 



energy is given by E, 






2.2n2 



IHK:K)p 



■ t E(°)(n*)-E(°)(n) 



rfn 



n* + 1 



U(n* - 1) - ix ix- Un* 



Therefore the energy spectrum E is given by the eigenval- 
ues of the Hamiltonian H Eq. ([SJ, i.e. 



2„,2 



with 



t(6)N 



1 + t- 



E = const + m 7 



n* n* + 1 

U(n* - 1) - ju p-Un* 



(12) 



(13) 



The phase transition between a Mott-insulator phase 
where 7 = and a supcrfluid phase where 7 > occurs 
when m = 0. Therefore the phase diagram at T = is 
given by 



t B (U,n,T = 0) = U 



(6) [n/U-n*][(n* 



1) 



m/c/] 



(^ 2 > 



M/t/ + 1 



(14) 

with |u/[/ e [n* — l,n*]. The difference with respect to 
the mean-field phase diagram for regular lattices is that t c , 
the critical hopping rate, depends on the second moment 
of the expected degree distribution, i.e. (6 2 ) = (k(k — 1)). 
Given the general expression for the expected degree dis- 
tribution of complex networks considered in this paper Eq. 
Q including the exponential cut-off £, the Mott-insulator 
phase disappears as £ — > 00 when A < 3 and remains finite 



instead if A > 3 and £ 



o. Therefore as -W- diverges 
i.e. as £ — > 00 while A < 3 wc have that the Mott insula- 
tor phase shrinks and finally disappears for large network 
sizes. Also it can be seen that the critical indices will de- 
viate from the mean-field values and they can be found 
by applying the heterogeneous mean-field techniques [6] 
developed for the classical phase transition. 

At finite temperature we cannot properly speak about 
a Mott insulator phase but wc have still a phase diagram 
between a normal phase and the supcrfluid phase. The 
local order parameter is given by the thermal average of 
the creation and annihilation operators, i.e. 



</v 



Traie- 0Hi 



(15) 



Using the same steps as in [29] we can prove that the crit- 
ical line for the Mott-insulator, supcrfluid phase is given 
by 



(0) 



r=0 e 



/3[/M--(£//2)r(r-l)] 



(11) t c (U,fX,p) (Q2)Y^ =Q Q r (lJ : n)e^-{U/2)r{T-l)] 



The ground state energy E\ \n) = E^°\n) with where 
£(°)(n*) = Oif/i < 0and£( o )(n*) = -pn* + ±Un*(n*-l) 
if p € (U(n* — 1), Un*) The second order correction to the 



MU,n) 



p + U 



(16) 



(17) 



(p-Ur)(U(r-l)-p) 

Therefore the phase diagram at finite temperature is also 
affected by the topology of the network and significantly 
changes when (0 2 ) diverges. 

Phase diagram of the Bose-Hubbard phase tran- 
sition on quenched complex networks. We now 

discuss the phase diagram of the Bose-Hubbard model on 
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Fig. 1: (Color online) Average order parameter for the super- 
fluid phase for scale-free networks with power-law exponent 
A = 2.2 and network sizes TV = 100 (top)TV = 1000 (mid- 
dle) and TV = 10,000 (bottom). As the network size increases 
the phase diagram changes monotonically as predicted by the 
mean-field treatment in the case A < 3. Therefore there is no 
Mott-insulator phase in the limit TV — > oo. 

quenched networks. In particular we focus on the base 
diagram at T — 0. The Hamiltonian we consider is the 
original Bose-Hubbard Hamiltonian defined on the adja- 
cency matrix r (Eq. ([I])). We solve this equation again 
in the mean-field approximation assuming Eq. (|6]) and 
•i/ji — (oj) = (a]). By solving self-consistently for ipf we 
found the boundary of the Mott-insulator phase where 
ipi = and the supcrfluid phase where ipi > OVz. The 
mean-field Hamiltonian H MF is then parametrized by the 
self-consistent parameters ipi and reads 



H 



MF 



E 



u 



i(m - l) - nm - 1 y^ Tjj (oj + al)ijjj 
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+'EE T ^*- 



(18) 



Following [3D] we consider the hopping term as a pertur- 
bation. At the first order of the perturbation theory we 
get that 

^ i = ^F{ i x,U)Y,n^ j (19) 

3 

where 



F{n,U) 



3 

n + u 



(20) 



[fj,-n*U][U{n*-l)-fj] 

where (i € (U(n* — 1), Un*). Therefore, if A is the max- 
imal eigenvalue of the adjacency matrix {ry}, the Mott- 
insulator phase with ipi = is stable as long as 



-F(jt,U)A<l. 



(21) 



Fig. 2: (Color online) Average order parameter for the su- 
perfluid phase for scale-free networks with power-law expo- 
nent A = 3.5 and network sizes TV = 100, TV — 1000 and 
TV = 10, 000. As the network size increases the phase diagram 
has slower finite size effects with respect to the case A < 3. 
Nevertheless the maximal eigenvalue increases with network 
size as A oc \/k max and therefore the Mott Insulator phase 
disappears in the thermodynamic limit also in this case. 



We observe that in random scale-free networks with de- 
gree distribution p(k) = Afk~ x the maximal eigenvalue A 
of the adjacency matrix diverges with a diverging value of 
the maximal degree of the network k max as A oc \Jk max 
48 51]. Therefore we find that also as long as the max- 
imal degree of the network diverges the Mott-insulator 
phase disappears in the large network limit, changing the 
phase diagram of the model significantly with respect to 
regular networks where the maximal degree of the net- 
work remains constant. We have checked these results 
by performing numerical integration of the mean-field cal- 
culations. We have studied the phase diagram of sin- 
gle quenched networks with scale-free degree distribution 
p(k) = Nk~ x and different values of the power-law expo- 
nent A to see how fast the convergence of the solution to 
the asymptotic phase diagram is. In the following we will 
show our finite-size scaling calculations and the resulting 
effective phase diagram of the Bose-Hubbard model within 
the mean-field approximation on the quenched network 
for different values of the number of nodes TV. In figure 
[T] we plot the effective phase diagram for network sizes 
N = 100, 1000, 10000 finding that for A = 2.2 < 3 the 
boundary of the Mott-insulating phase decreases with the 
network size. On the other hand for a typical network with 
A = 3.5 > 3 (sec figure |2J) the phase diagram has slower 
finite size dependencies. This shows that the annealed ap- 
proximation for the Bose-Hubbard model on scale- free net- 
works strongly differs from the quenched phase diagram of 
the model for 7 > 3. In particular we have that the Mott 
insulator phase transition on annealed scale- free networks, 
with diverging second moment of the degree distribution, 
vanishes in the thermodynamic limit while in quenched 
networks it is sufficient that the most connected node has 
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Fig. 3: (Color online) Left panel: the first 3 generations of the 
Apollonian graph. The nodes and links added to construct the 
fst (red), 2nd (green) and 3rd (purple) generations are shown. 
Right panel: the 5th generation Apollonian network. 



a diverging connectivity in the thermodynamic limit to 
destroy the Mott insulator phase. 

Finally we have studied the phase diagram as predicted 
by the mean-field approximation, on Apollonian networks 
[20] . Apollonian networks are an example of scale- free 
networks which are embedded in a two dimensional space. 
They are constructed by considering the classical 2D Apol- 
lonian packing model in which the space between three 
tangent circles placed on the vertices of an equilateral 
triangle is filled by a maximal circle. The space-filling 
procedure is repeated for every space bounded by three 
of the previously drawn tangent circles. The correspond- 
ing Apollonian network is constructed by connecting the 
centers of all the touching circles (Fig. [3] (left)). The 
resulting network is scale-free with power-law degree dis- 
tribution p(k) = Afk~ x and A = 1 + ln(3)/7n(2) ~ 2.585. 
Also these networks are known to have diverging maximal 
eigenvalue A of their adjacency matrix [21]. Therefore 
we expect that also in these networks the Mott-insulator 
phase of the Bosc-Hubbard model should disappear in the 
large network limit. We have used the iterative scheme in- 
troduced in [21] to generate 5th (Fig. [3] (right)), 7th and 
9th generation Apollonian networks, where the number of 
nodes in the nth generation is N(n) = (3"+5)/2. In Fig. |4] 
we demonstrate the finite-size effects on the phase diagram 
found in the mean-field approximation. We note here that 
in experimental realizations of the Bose-Hubbard model 
on Apollonian networks it might be relevant to include a 
hopping term that depends also on the distance between 
the nodes. The dynamics on this weighted network might 
further modify the phase diagram of the process. In this 
paper we have chosen to study the simple case in which 
the hopping doesn't depend on the distance, leaving the 
characterization of the dynamics on the weighted network 
for future investigations. 

Conclusions. — In conclusion in this paper we have 
shown that the scale- free network topology of the underly- 
ing network strongly affects the phase diagram of the Bose- 
Hubbard model. By performing mean-field calculations on 
annealed networks we have shown that the Mott-insulator 
phase disappears in the large network limit as long as 
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Fig. 4: (Color online) Effective phase diagram of the Bose- 
Hubbard model on Apollonian networks of the 5th, 7th and 
9th generation, with network sizes N = f24, N — f096 and 
N = 9844, respectively. As the network size diverges the Mott- 
insulator phase is reduced and it disappears in the limit of an 
infinite network. 



the power-law exponent A of the degree distribution is 
A < 3 and the exponential cutoff £ of the distribution di- 
verges. We have performed mean-field calculations in an- 
nealed networks finding in this approximation the phase 
diagram of the model both at T = and at finite temper- 
ature. Moreover the analytical and numerical solutions of 
the Bosc-Hubbard model on quenched networks show that 
this argument must be corrected in the quenched case and 
that it is sufficient that the maximal eigenvalue diverges in 
order to change the phase diagram of the model. Finally 
we have considered the Bose-Hubbard model on Apollo- 
nian networks that are an example of scale-free networks 
embedded in a two dimensional space. In short, this work 
offers a new perspective on the characterization of quan- 
tum critical phenomena in annealed and quenched com- 
plex networks and shows that the second moment of the 
degree distribution (k 2 ) and the maximal eigenvalue of 
the adjacency matrix play a crucial role in determining 
the phase diagram of the Bose-Hubbard model. 
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